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Ph' Abstract 

Meson-charmonium dissociation reactions governed by the quark interchange are 



studied with temperature-dependent quark potentials. Quark-antiquark relative- 



's- , motion wave functions and masses of charmonia and charmed mesons are deter- 

o 

mined by the central spin-independent part of the potentials or by the central spin- 
independent part and a smeared spin-spin interaction. The prominent temperature 

\6 . 

dependence of the masses is found. Based on the potentials, the wave functions, 



and the meson masses, we obtain temperature-dependent cross sections for the fif- 
teen dissociation reactions: trJ/tp —> D*D or DD*, nJ/ip —> D*D* , Trip' — > D*D 



or DD*, nip' -> D*D*, irxc -> D*D or DD*, irxc -»• D*D*, pj/ip -»• DD, 
pj/ip -»• D*D or DD*, pj/ip -> D*D* , pip' -»• DD, pip' -»• D*D or DD*, 
pip' -»• /)x c -»• 5 A px c D*D or DD*, and pxc -> D*D* . The nu- 

merical cross sections are parametrized for future applications in hadronic matter. 
The particular temperature dependence of the J/ip bound state leads to unusual 
behavior of the cross sections for endothermic J/ip dissociation reactions. The quan- 
tum numbers of ip' and Xc can not make their difference in mass in the temperature 
region 0.6T C < T < T c , but can make the ip' dissociation different from the Xc 
dissociation. 

PACS: 25.75.Nq, 12.39.Jh, 13.75.Lb 
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1. Introduction 

From the year when Matsui and Satz [lj originally suggested the suppressed J /if) 
production as a signature for the formation of a quark-gluon plasma in high-energy heavy- 
ion collisions to the last year's quark matter conference where preliminary J/if> data in Pb- 
Pb collisions at the Large Hadron Collider were reported by the ALICE Collaboration [2] 
and the CMS Collaboration [3], important measurements have been done. Essential 
theoretical progress on J/if> has also been made. One of the fundamental issues is the 
dissociation of charmonia in hadronic matter [I]. To identify J /if} as a probe of the quark- 
gluon plasma in a definite way, hadron-charmonium dissociation processes must be well 
understood. Calculations of dissociation cross sections are thus an important aspect in 
studying J/if> physics. The dissociation is described by the meson or quark degree of 
freedom. Corresponding to the two degrees of freedom, different scattering mechanisms 
can be assumed, and different results on the dissociation cross sections have been reported 
in the literature. 

There are mainly three approaches to the meson- J /if} dissociation problem. In the 
short-distance approach the parton model of light hadrons, the constituent quark model of 
J /ip, and the ghion- J/if) dissociation cross section of Peskin and Bhanot [5] are employed in 
Refs. P,[7j to investigate nJ/if) and NJ/if> dissociation. The portion of hard gluons inside 
the nucleon or pion at low energies is not large enough to induce mb-scale cross sections. 
In the meson-exchange approach J /if} dissociation has been studied with effective meson 
Lagrangians. Since the discovery of J /if) reveals the charm quark, J /if} is first described by 
QCD. The meson field description of J /if) and charmed mesons began to appear from the 
work of Matinyan and Miiller [8] who first discussed the t-channel exchange of D meson 
in inelastic ivJ/ip and pJ/ip scattering. They got mb-scale cross sections for irJ/ip — y 
D*D + DD* and pJ/if) — > DD at low energies. A similar scale of cross sections has 
been obtained in other meson Lagrangians with different symmetries and modified vertex 
functions in Feynman diagrams to include the effect of finite meson form factors [9HT4]. In 
the quark-interchange approach J /if) dissociation has been studied in nonrelativistic quark 



3 



potential models pTSHTH] . The assumption that color- independent confining interaction 
acts only between a quark and an antiquark in Ref. [T5] yields that the J /if) absorption 
cross section corresponding to nJ /if) — > D*D + DD* + D*D* has a peak value of about 
7 mb at the kinetic energy = \fs — (m„ + mj/^f) ~ 0.8 GeV. However, by including 
color generators in the linear confining potential and allowing the interaction to connect 
any two constituents (quarks and antiquarks), Wong, Swanson, and Barnes [16] gave a 
peak value of only ~ 1 mb at the same kinetic energy. Usefully, the dissociation cross 
sections of ground state, orbitally and radially excited charmonia in collisions with 7r and 
p mesons have been presented in Ref. [18] where all parameters in the color Coulomb, spin- 
spin hyperfine and linear confining interactions are determined by fits to the experimental 
meson spectrum. 

All the above works only involve charmonium dissociation reactions in vacuum, and 
have been used to evaluate J /if) suppression in nucleus-nucleus collisions, but we know 
that hadron masses, the quark potential and so on are affected by medium and the QCD 
phase transition is the most striking medium effect, so we must study medium effects 
on the charmonium dissociation. In this work we calculate dissociation cross sections of 
J /if), if)' i and Xc i n collisions with tc and p at various temperatures. The quark-interchange 
mechanism [19], the Born approximation, and a temperature-dependent quark potential 
[20] are ingredients in establishing cross section formulas. 

This paper is organized as follows. In the next section we use the temperature- 
dependent potential [20] in the Schrodinger equation to obtain temperature-dependent 
masses of charmonia and charmed mesons. In Section 3 we give formulas for charmonium 
dissociation cross sections. Numerical results for unpolarized cross sections of J /if), if)', and 
Xc in collisions with 7r and p at six temperatures are shown in Section 4 and relevant dis- 
cussions are given. In Section 5 the spin-spin interaction arising from one-gluon exchange 
is smeared. Subsequent results are shown and numerical cross sections are parametrized. 
In Section 6 we present a procedure on how to obtain unpolarized cross sections at any 
temperature in the region 0.65 < T/T c < 1 where T c is the critical temperature. Finally, 
we summarize the present work in Section 7. 
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2. Masses of charmonia and charmed mesons 



As shown in Refs. [2ITI24"] . the quark potential containing flavor-independent confine- 
ment, a Coulomb term and hyperfine interactions can consistently reproduce masses from 
light to heavy hadrons, and the flavor-independent assumption of confinement is thus 
reasonable. Such confinement in hadronic matter can be estimated by the lattice calcula- 
tions [25J and depends on temperature. At large distances the confinement manifests itself 
by a plateau that lowers with increasing temperature. In Ref. [20] we used the confinement 
at large distances and the short-distance potential originating from one-gluon exchange 
and loop corrections in perturbative QCD [26] to construct a central spin-independent, 
flavor-independent but temperature-dependent potential 



Vsi(f) = • D 

bi\ j 2 2 4 



1.3 



T N 4 



tanh(Ar) + ^.^!^exp(-£r). (1) 



T 

Only here D = 0.7 GeV, T c = 0.175 GeV, A = 1.5[0.75 + 0.25( J) 10 ] 6 GeV, E = 0.6 GeV, 
and A = y^3b /lQn 2 a' in which a' = 1.04 GeV" 2 is the Regge slope and b = 11 — |iV/ 
with the quark flavor number Nf = 4. A a are the Gell-Mann matrices for the color 
generators of constituent a. The dimensionless function v (x) is an integration over the 
absolute value of gluon momentum Q 

Ab [°°dQ, _,- 2 



v{x) = — / -%(p(Q*) - -=r) sin(^), (2) 
^ Jo Q Q 2 A 

where K = 3/167r 2 a'. The subtraction of K/Q 2 from the physical running coupling 

constant p{Q 2 ) leaves only the contribution of one-gluon exchange plus perturbative one- 

and two-loop corrections. The factor exp(— Er) is a medium modification factor to the 

potential of one-gluon exchange plus perturbative one- and two-loop corrections. The 

temperature dependence is completely negligible at very short distances and obvious at 

intermediate and large distances. The potential well fits the lattice gauge results at 

T/T c > 0.55 [25], but fails to give a long-ranged linear confining potential at T = 0. 

Given the charm quark mass m c = 1.51 GeV, the Schrodinger equation with the 

potential given in Eq. (1) is solved to obtain masses and wave functions of J/if>, if>', and 

Xc- Here the \c mass corresponds to the center of gravity of Xco, X&, an d Xc2 [27J. At 



T = the masses of J/if), f, and Xc are 3.10505 GeV, 3.67679 GeV, and 3.51138 GeV, 
compared to the experimental values 3.096916 GeV, 3.68609 GeV, and 3.5253 GeV [25] . 
respectively. The temperature dependence of the charmonium masses is shown by the 
solid, dashed and dotted curves in Fig. 1. Because J/ip has a small radius, the J/ip mass 
changes more slowly than the ip' and \ c masses. Even though ip' and \ c have different 
quantum numbers, both become degenerate in mass in the temperature region in Fig. 1. 
Furthermore, at temperatures very close to the critical temperature J/ijj joins and Xc 
to form a triplet in mass. This is a medium effect on charmonia! 

Given m c = 1.51 GeV, the up and down quark masses m u = =0.32 GeV, and 
the strange quark mass m s =0.5 GeV, the Schrodinger equation with the central spin- 
independent potential offers the same quark-antiquark relative-motion wave functions of 
it and p (D and D*, D s and D*) and the spin-averaged mass of ir and p (D and D*, D s and 
D*). The spin-averaged mass of a spin-0 meson and a spin-1 meson with the same isospin 
is one-fourth of the spin-0 meson mass plus three-fourths of the spin-1 meson mass. The 
quark-antiquark relative- motion wave functions of tt and p (D and D*, D s and D*) are 
used to calculate the mass splitting of tt and p (D and D*, D s and D*) with the spin-spin 
interaction that arises from one-gluon exchange plus one- and two-loop corrections [29] 




where s a and m a (nib) are the spin and mass of the constituent a (b), respectively. 
From the mass splitting and the spin-averaged mass we get the spin-0 meson mass and 
the spin-1 meson mass: the former is the spin-averaged mass minus three-fourths of the 
mass splitting; the latter is the spin-averaged mass plus one-fourth of the mass splitting. 
At T = the masses of D, D*, D s , and D* are 1.80666 GeV, 2.09552 GeV, 1.85228 
GeV, and 2.18695 GeV, compared to the measured values 1.86722 GeV, 2.00861 GeV, 
1.96847 GeV, and 2.1123 GeV, respectively. The temperature dependence of the D, D*, 
D s and D* masses is plotted in Fig. 1 as two long dashed curves and two dot-dashed 
curves. The D and D s masses keep almost unchanged from T = 0.6T C to 0.85T C and 0.9T C , 
respectively, and apparently fall off in the other temperature regions. The D* and D* 




ss 



Xa Xb 167T 2 

2 ' 2 25 



S 3 (r) 



s, a ■ s b | A a Afc 47r 1 d 2 v(Xr) s a ■ s b 
m a mb 2 2 25 r dr 2 m a mb ' 



(3) 
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masses decrease slowly from T = 0.6T C , and become apparently falling off from T = 0.8T C 
and T = 0.85T C , respectively. The medium effect on charmed mesons is obvious only in 
the region where the masses apparently fall off. The temperature dependence of ir and 
p masses was shown in Fig. 2 of Ref. [20]. With increasing temperature the ir mass 
decreases slowly for 0.6T C < T < 0.78T C and rapidly for T > 0.78T C while the p mass 
decreases rapidly for 0.6T C <T <T C . From T = 0.6T C to 0.99T C the masses of ir, p, J /if), 
ip', Xc, D, D*, D s , and D* are reduced by 100%, 99%, 7%, 16%, 16%, 20%, 27%, 12%, and 
23%, respectively. Therefore, the medium effect on the two light mesons is more obvious. 
Either D and D* or D s and D* become a doublet in mass at T — > T c , but the four mesons 
do not become degenerate unlike J/ip, ip' and Xc- 

The meson masses in units of GeV in the region 0.6 < T/T c < 1 are parametrized as 

0.03 



m 



J/4> 



3.07 



1 



T 



m$> = 3.48 



m 



Xc 



3.42 



1.01T t 

T 
2.19T C 

T 
1.89T r 



m D = 1.795 



T 



m D * 



2.02 



1.94 



2.133 



1 - 



1.16T C 

T 
1.42T r 



T 



1.02T C 

T 
1.29T r 



3.76 



4.63 



5.65 



9.67 



5.38 



3.3 



7.74 



6.96 



0.92 



2.18 



0.08 



6.28 



1.29 



(4) 

(5) 
(6) 
(7) 
(8) 
(9) 
(10) 



3. Cross section formulas 



In nonrelativistic dynamics for the quark-interchange process qq + cc — )• qc + cq, the 
center-of-mass motion of qq and cc (i.e. qc and cq) is separated off. This guarantees 



that cross sections are calculated in a way independent of the center of mass. We then 
choose the center-of-mass frame where the cross section can be easily formulated [3"0] . 
We denote the mass and the four-momentum of meson i (i = qq, cc, qc, cq) by irti and 
Pi = (Ei, Pi), respectively The Mandelstam variables are s = (E qq + E c5 ) 2 — (P qq + P C c) 2 
and t = (E qq — E q5 ) 2 — (P qq — P q5 ) 2 . The cross section for the meson- charmonium scattering 
qq + cc — > qc + cq in the center-of-mass frame is 

a(S,m s ,V~s,T) = 1 Iffigi de\M R (s,t)\ 2 sine, (11) 
32tts \P(y/s)\ Jo 

where S is the total spin of either the two incoming mesons or the two outgoing mesons, 
ms is the magnetic projection quantum number of S, is the transition amplitude, 
and 9 is the angle between qq momentum P and qc momentum P'. P and P' are related 
to the Mandelstam variable s by 

P\^) = ^{[s- K- + m%)] 2 - 4m^-m c 2 c-} , (12) 

= Ts { t s ~ ^ + ~ Am ^} • (13) 

The interchange of quarks brings about two scattering forms, the prior form and the 
post form. The two forms may lead to different values of the transition amplitude .Ms 
(and the subsequent cross section), which is the so-called post-prior discrepancy [3TH33] . 
Scattering in the prior form means that gluon exchange takes place prior to the quark 
interchange, and the corresponding transition amplitude is 

Mf ov = A^/E qq E c - c E q - c E cq < ^ q5 \ < Vcs| (V g - C + V cq + V qc + V q5 )\^ qq > |^ cc - >, 

(14) 

while scattering in the post form means that the quark interchange is followed by gluon 
exchange, and the corresponding transition amplitude is 

Mfi° St = ^E qq E c ,E q - c E cq < ^ q5 \ < ip cq \{V qq + V c5 + V qc + Vqc) \^ qq - > |^ cc - >, 

(15) 

where ipq 5 is the product of color, spin, flavor and momentum-space wave functions of the 
relative motion of q and c and satisfies J j5^^qc(Pqc)4>qc(Pqc) = 1 where p q c is the relative 
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momentum of q and c, and similarly if) cq -, ip qq - and ip c5 . A relative momentum depends on 
a linear combination of P and P'. The momentum-space wave functions are the Fourier 
transform of the coordinate-space wave functions which are solutions of the Schrodinger 
equation. In the transition amplitude we use the Fourier transform of the sum of the 
central spin-independent potential and the spin-spin interaction: 



1.3 



\Z&tX\ Sn r ^ rsin(Qr) 



(2n) d 5 d (Q) - — / dr 



Q J exp(2Ar) + 1 



2 2 J (E 2 + Q 2 + q 2 ) 2 - AQ 2 q 2 

X a X b 16vr 2 s a ■ s b X a X b 16tt 2 A f°° d 2 v (x) . ( Q \ s a ■ s b 
+ tt ■ — ^ / dx — — — sin —x 



2 2 25 m a m b 2 2 25Q J dx 2 \X J m a m b 

(16) 

Let cr pnor and a post be the cross sections for scattering in the prior form and in the post 
form, respectively, and they are given by Eq. (11). The unpolarized cross section for 
qq + cc — >■ qc + cq is 

1 



a unpol (y/ r s, T) 



(2S qq +l)(2S c - c + l) 

x J2(2S + 1) (TPri ° r( ^ ms ' v/i ' T) + ^ (g ' m5 ' ^ T) 
s 2 

(17) 

where S qq - and S'cg are the spins of qq and cc, respectively. 

4. Numerical cross sections and discussions 

Quark-antiquark relative-motion wave functions of mesons are given by the Schrodinger 
equation with the central spin-independent potential V si . With the up (down) quark mass 
0.32GeV and the experimental tt mass, the experimental data of 5-wave 1 = 2 elastic 
phase shifts for scattering in vacuum [341437] for < y/s < 2.4 GeV are reproduced 
from the potential in Eq. (16) at T = (20]. In the estimates of charmonium dissociation 
cross sections at T = 0, the experimental masses of pion, rho, charmonia, and charmed 
mesons are employed. For 0.6 < T/T c < 1 we use the temperature-dependent meson 
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masses shown in Fig. 1 and the ir and p masses in Fig. 2 of Ref. [20]. Temperature de- 
pendence of charmonium dissociation cross sections shown in the next subsection seems 
to be complicated, but is understandable. 

From T = 0.6T C to 0.99T C the masses of tc, p, J/tp, ip', Xc D, and D* are reduced by 
0.40955 GeV, 0.60936 GeV, 0.21908 GeV, 0.56022 GeV, 0.53639 GeV, 0.35318 GeV, and 
0.53845 GeV, respectively. The reduced amounts of ip' , Xc an d D* masses are between 
the ones of 7r and p masses. Hence, the difference m q5 + m cq - — m q q — m c5 may be larger or 
smaller than zero. Whether a reaction is endothermic or not depends on temperature. A 
reaction may be exothermic below a temperature and endothermic above the temperature. 

4.1. Numerical cross sections 

The temperature-dependent potential given in Eq. (16), the wave functions, and the 
meson masses in the transition amplitude make both the meson-charmonium dissociation 
cross sections and relevant threshold energies dependent on temperature. In Figs. 2-16 
we plot cross sections for the following fifteen meson-charmonium dissociation reactions: 
TrJ/ip ->- D*D or DD*, nJ/ip ->• D*D*, nip' -> D*D or DD*, nip' ->- D*D*, 7i Xc -> D*D 
or DD*, nxc -»■ D*D*, pj/^j DD, pJ/4> -»■ D*D or DD*, pj/tp D*D*, pip' DD, 
pip' -»■ D*D or DD*, pip' D*D*, p Xc -»■ DD, PXc -»> D*D or DD*, and px c -»■ D*D*. 
Here D stands for D+ or D°, D for D" or D°, D* for D*+ or D*°, and D* for D*~ or D*°. 
For the reaction that produces D*D or DD*, the cross section for the production of DD* 
equals the one of D*D. No matter what the temperature is, the reactions, nJ/ip — > D*D, 
nJ/ip ->■ D*D*, Trip' -> D*D*, nx c -> ^D*^*, pJ/V> ->■ an d -> ^*^*, are 

endothermic. No reactions are all exothermic from T = 0.6T C to T c . To display in- 
medium modification of dissociation, we plot the cross sections at T = for reference 
even though they differ from the results obtained in Ref. [18] from the potential that 
includes the color Coulomb, spin-spin hyperfine and linear confinement terms. 

4.2. Endothermic reactions of Trip' , ttXc pip' and px c 
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The cross sections for endothermic nip' and nx c reactions in Figs. 4-7 show that 
peak cross sections increase from T/T c = to 0.75 and decrease from T/T c = 0.75 to 
0.95 except the increasing peak cross section from T/T c = 0.9 to 0.95 in Fig. 5 for 
nip' — y D*D*. While temperature increases from zero, the large-distance potential given 
by the first term in Eq. (1) becomes smaller and smaller, and the Schrodinger equation 
with the potential shown by Eq. (1) produces increasing meson radius. The increase of 
initial-meson radii leads to peak cross sections increasing from T/T c = to 0.75. With 
increasing temperature, it becomes more difficult to combine the interchanged quarks and 
the antiquarks to form mesons qc and cq due to weakening confinement. This leads to 
peak cross sections decreasing from T/T c = 0.75 to 0.95. This also explains the decrease 
of the peak cross section from T/T c = 0.85 to 0.9 in Fig. 13 for pip' — y D*D*, from 
T/T c = 0.9 to 0.95 in Fig. 11 for pip' ->■ DD, Fig. 12 for pip' -> D*D, Fig. 14 for 
PXc ->■ DD and Fig. 15 for px c ->• D*D, or from T/T c = 0.85 to 0.95 in Fig. 16 for 
p Xc -> D*D*. 

We have already known that ip' and Xc have very similar temperature dependence in 
mass. The nip' (pip') and nx c (PXc) reactions are almost identical in threshold energy, 
but are different in cross section even though the final mesons in the nip' (pip') reactions 
are the same as those in the nx c (PXc) reactions. The difference is clearly attributed 
to the different quantum numbers of ip' and Xc- Particularly, when T/T c goes from 0.9 
to 0.95, the increase of the peak cross section of nip' — y D*D* in Fig. 5 (pip' — y D*D* 
in Fig. 13) is against the decrease of the peak cross section of nx c — > D*D* in Fig. 7 
(PXc — > D*D* in Fig. 16). This relates to the node in the ip' wave function. The node 
leads to cancellation between the negative wave function on the left of the node and the 
positive wave function on the right of the node in the integration involved in the transition 
amplitude. While the cancellation at T = 0.95T C is less than at T = 0.9T C , the peak cross 
sections of nip' ->■ D*D* and pip' ->■ D*D* rise from T/T c = 0.9 to 0.95. 

4.3. Exothermic reactions of nip', nx c , pip' ' , PXc and pJ /ip 

Since cross sections for exothermic reactions at threshold energies are infinite, we start 
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calculating the cross sections at y/s = m qq + m c5 + A^/s with Ay's = 10~ 4 GeV. At the 
energies the cross sections correspond to the curve tops, for example, the top of the curve 
for T/T c = 0.65 displayed in Fig. 4, and from Eq. (12) we have 

p2~2Av^/(— + — ), (18) 

which indicates that P 2 decreases with decreasing meson mass. In Section 2 we have ob- 
tained that the meson masses decrease with increasing temperature. Hence, \P\ decreases 

— * o 

while temperature increases. At T = 0.6T C , \P\ is estimated at about GeV. The 

— * 

influence of P on the relative momenta p a b is thus negligible. Finally, p a b is governed by 
P' that obeys 

m qq + m c5 + Ays = \l m 2 q - c + P? + \J ' m 2 cq + P' 2 [ (19) 
Since Ay's is very small, we have 

P' w 2(m qq + m cd - m q c - m C q) / { 1 ). (20) 

With increasing temperature the value of m qq + m c5 — m q5 — m cq decreases in the region 
where reactions are exothermic, and becomes zero at a temperature; meanwhile, \P'\ 
decreases so that the absolute value of the relative-motion part of ip a b(Pab) gets larger 



and larger, and the factor E qq E ce E qe E cq in the transition amplitude decreases. Simple 
calculations show that the factor generally decreases with increasing temperature in 
the region where reactions are exothermic. The product of the increasing absolute value 
of the relative- motion part of ip a b{p a b) and the decreasing ^j^E qq E c5 E q cE cq leads to that 
the cross section at \/s = m qq + m c5 + Ay/s first increases and then decreases as shown 
in Figs. 6, 8, and 16, or first decreases and then increases in Figs. 4 and 11 through 15. 

4.4. Endothermic reactions of nJ/ip and pJ/ip 

As shown in Figs. 2-3 and 8-10, the peak cross sections of the endothermic nJ/ip and 
pJ/tp reactions decrease from T = to 0.85T C and increase from T = 0.85T C to 0.95T C . 
The behavior of the peak cross sections is different from that of the endothermic reaction 
nxc D*D*, which is an increase from T = to 0.75T C and a decrease from T = 0.75T C 
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to 0.95T C . The behavior must come from a special point, and as demonstrated below, the 
special point is just the temperature dependence of the mass and radius of J/ip. 

First, we discuss cross sections for nJ/ip and pJ/ip dissociation in the region < 
T/T c < 0.85. In Section 2 we have obtained the slow decrease of ip', Xc and D* masses 
and the very slow decrease of J/ip and D masses when T/T c increases from to 0.85. It 
is shown in Fig. 2 of Ref. [20] that the ir mass decreases slowly for 0.6T C < T < 0.78T C 
and the p mass does rapidly. Then, keeps roughly unchanged or decreases. Since the 
J/ip radius changes very slowly in comparison to the increase of the it or p radius, only 
one increasing initial-meson radius causes the increase of |A^fi| 2 , and the increase can not 
overcome the amount reduced by the weakening confinement with increasing temperature, 
i.e. |.Mfi| 2 decreases. Finally, the product of and |.Mfi| 2 leads to the decrease of the 
peak cross sections of the endothermic irJ/ip and pJ/ip reactions when T/T c increases 
from to 0.85. 

Next, we discuss cross sections for nJ/ip and pJ/ip dissociation in the region 0.85 < 
T/T c < 1. The factor increases by about 50% from T = 0.85T C to 0.9T C and by about 
100% from T = 0.9T C to 0.95T C . When temperature increases from T = 0.85T C to 0.95T C , 
not only the it and p radii increase, but also the J/ip radius apparently does. Now the 
increase of |.Mfi| 2 caused by the increasing radii of the two initial mesons overcomes the 
amount reduced by the weakening confinement with increasing temperature. Finally, the 
peak cross sections of the endothermic nJ/ip and pJ/ip reactions increase from T/T c = 
0.85 to 0.95. 

4.5. Comparison ofn+ charmonium reactions and p+ charmonium reactions 

We make a comparison of peak cross sections of endothermic it + charmonium reac- 
tions with ones of endothermic p + charmonium reactions. We find that the peak cross 
section of nJ/ip — > D*D (tcx c — > D*D) is larger than that of pJ/ip -» D*D (px c — > D*D) 
at the same temperature while the peak cross section of 7r + charmonium — )■ D* + D* 
(7T0' — > D*D) is smaller than the one of p + charmonium — > D* + D* {pip 1 — > D*D). 
Exothermic tt + charmonium reactions and exothermic p+ charmonium reactions are 
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compared about their cross sections at y/s = m qq + m cg + 10" 4 GeV (the largest cross 
section shown in each curve). It is shown in Figs. 4, 6, 12 and 15 that the cross sec- 
tion of 7T+ charmonium —)>/)* + D at T/T c = 0.65 or 0.75 is larger than the one of p+ 
charmonium —)•/)* + D at the same temperature. It is meaningless to compare endother- 
mic 7r + charmonium reactions (e.g. m/j' — > D*D at T/T c = 0.85) with exothermic p + 
charmonium reactions (e.g. pip' — > D*D at the same temperature). 

Regarding endothermic reactions we examine P and P' at y/s = m q5 + m cq + d at 
which the peak cross sections are given and do <C m q5 + m cq . Eq. (13) yields 

P' 2 w 2d /(— + — ). (21) 

From an endothermic 7r + charmonium reaction to an endothermic p + charmonium 
reaction at the same temperature, m q5 and m cq do not change, and do decreases or does 
not change. Accordingly, P' is reduced or unchanged, is related to initial-meson 
energies by 

m q - c + m cq + d = sJm 2 qq + P 2 + sJm^ + P 2 , (22) 

which let P 2 decrease from the endothermic it + charmonium reaction to the endothermic 
p + charmonium reaction. As a result the absolute value of the relative-motion part of 
ipabiPab) in the transition amplitude increases. 

In the transition amplitude ip a b is the product of color, spin, flavor and momentum- 
space wave functions of the relative motion of constituents a and b. The transition am- 
plitude involves the spin factor: the matrix elements of s a ■ s b of the spin-spin interaction 
and the overlap of the spin wave function of final mesons and the one of initial mesons. 
The difference of the transition amplitudes for tt + charmonium — > D* + D and for p + 
charmonium —)>/)* + D is caused by not only the difference of the ir and p masses (the 
subsequent changes of yj E qq E c5 E q5 E cq , ip a b(Pab), and \P'\/s\P\) but also the difference of 
the spin factors for the two reactions. Particularly, n and p become almost degenerate 
in mass at T/T c > 0.9. Then, only the spin factors made the peak cross section of n + 
charmonium — >■ D* + D larger or smaller than the one of p + charmonium — >■ D* + D. 
In case the final mesons are D* and D*, the total spin of ir and charmonium is 1 that 
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differs from the total spins 0, 1 and 2 of p and charmonium. The spin factors plus the 
additional contribution of the channels of the total spins and 2 of p and charmonium 
make the peak cross section of tt + charmonium — > D* + D* smaller than the one of p 
+ charmonium —?►/)* + D* at T/T c = 0.9 and 0.95. The difference between the cross 
sections at y/s = m q q + m c5 + 10~ 4 GeV for exothermic tt + charmonium — > D* -\- D and 
for exothermic p + charmonium — )■ D* + D can be understood similarly with Eqs. (18) 
and (19). 

5. Results pertinent to smearing in the spin-spin interaction 

We have solved the Schrodinger equation with the central spin-independent potential 
V^i in Section 2. If we use the spin-spin interaction given in Eq. (3) in the Schrodinger 
equation, the delta function in the first term of the V ss expression cannot be correctly 
dealt with. However, in the present section we take smearing in the spin-spin interaction 
by the regularization 5 3 (f) — > exp(— d 2 r 2 \ and the smeared spin-spin interaction 

V - A b 167r 2 d 3 2 2 s a ■ s b X a \ b 47T 1 d 2 v(Xr) s a ■ s b 

ss ~ 2 ' 2 25 7r 3 / 26XP ^ T) m a m b 2 ' 2 25 r dr 2 m a m b ^ ' 

can be used in the Schrodinger equation. Such smearing actually includes relativistic 
effects [T71ES1EB]- The quantity d is related to quark masses by 

2 



a = cTq 



1 1 ( Am a m b 



2rn a m b 
m a + m b 



2 2\{m a + m b )\ 

where do = 0.15 GeV and o\ = 0.705. By solving the Schrodinger equation with the 
central spin-independent potential and the smeared spin-spin interaction, we obtain meson 
masses and quark-antiquark relative-motion wave functions that differ by different mesons. 
The p wave function is near that obtained by the Schrodinger equation with only V S i, but 
the it wave function is not. Since the smeared spin-spin interaction is used, the pion 
becomes a tight bound state of a quark and an antiquark. Then, the tt radius increases 
slowly from 0.6 to 0.9 of T/T c and quickly from 0.9 to 1 of T/T c . 

The mass splittings at T = are m p —m n = 0.6294 GeV and m^-»— itlk = 0.39865 GeV 
which are closer to the experimental data 0.6304 GeV and 0.3963 GeV than the results 
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(0.5989 GeV and 0.3833 GeV) [20] of the Schrodinger equation with only the central spin- 
independent potential. At T=0 we also get 3.13509 GeV, 3.69248 GeV, 3.50578 GeV, 
1.90578 GeV, 2.05274 GeV, 1.9614 GeV, and 2.13804 GeV as the masses of J/ip, ip', Xc 
D, D*, D s , and D*, respectively, compared to the measured values 3.096916 GeV, 3.68609 
GeV, 3.5253 GeV, 1.86722 GeV, 2.00861 GeV, 1.96847 GeV, and 2.1123 GeV [28]. For 
0.6 < T/T c < 1 temperature dependence of the p, K*, J/ip, ip', Xc D, D* , D s and D* 
masses obtained with the potential V B i + V^ s is very close to that obtained with only V s {. 
Therefore, the parametrizations given in Eqs. (4)-(10) are also valid in the present section, 
and the p and K* masses in units of GeV in the region 0.6 < T/T c < 1 are parametrized 
as 



nip = 0.73 



1 - 



T 



m K , 



0.84 



0.992T C 
T 



3.67 



0.989 



1.057; 



4.16 



(24) 



(25) 



However, the tt, K and rj masses obtained with V^; + V ss are smaller than those obtained 
with only V^. They are plotted in Fig. 17 for 0.6 < T/T c < 1 and are parametrized as 

0.51 



0.24 



T 



m K = 0.46 



0.55 



1 - 



1 



0.97T C 
T 

imt c 

T 

Tmt, 



3.81 



8.58 



3.11 



0.88 



0.29 



(26) 

(27) 
(28) 



in units of GeV. 

Corresponding to the tt mass, the reactions nJ/ip — > D*D, nJ/ip — > D*D*, Trip' — > 
D*D, nip' — > D*D*, nxc — > D*D, and iiXc — > D*D* are all endothermic. Cross sec- 
tions for the reactions are calculated with the experimental meson masses at T = 0, 
the temperature-dependent meson masses for 0.6 < T/T c < 1, the quark-antiquark 
relative-motion wave functions and the Fourier transform of the sum of the central spin- 
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independent potential and the smeared spin-spin interaction 



V ab (Q 



_A a A;, 3 

' ~2~i 

2 2 

A a A fe 167r 2 



2 25 



exp 



A a Xb 167r 2 A 
+ T ' ~2~ 25Q Jn 



7T 



(2vr)^(Q) - ^ 



r sin(Qr) 



Q Jo exp(2Ar) + 1 



„2 



(£ 2 + Q 2 + g 2 ) 2 - 4Q 2 g 2 

Q 2 \ S • S b 



4<i 2 / m a mb 

' d 2 v (x) . ( Q \ s a - s b 

ax — — — sin —x . 

ax 2 \ A / m a mb 



(29) 



We plot the cross sections for the reactions in Figs. 18-23. The cross sections differ from 
those shown in Figs. 2-7. For example, the reactions mp' — >■ D*D and -n\c D*D at 
T/T c = 0.65 and 0.75 are exothermic in Figs. 4 and 6 but endothermic in Figs. 20 and 
22. One feature for each reaction is that the peak cross section at T/T c = 0.95 is always 
larger than the ones at the other five temperatures T/T c = 0, 0.65, 0.75, 0.85, 0.9. This is 
caused by the 7r and J /if} radii which increase faster from 0.9 to 1 of T/T c than from 0.6 
to 0.9 of T/T c . The numerical cross sections for endothermic reactions are parametrized 
as 



a 



unpol 



(>/S,T) 



+a 2 



exp 
exp 



C 2 ( 1 



6i 



(30) 



where y/so is the threshold energy, and a 1 , 6 1; ci, a 2 , 6 2 and c 2 are parameters. Determi- 
nation of parameter values needs time-consuming computations, and the values are listed 
in Tables 1-3. 

Since the p, J /if), if)', Xc D and D* masses determined by the Schrodinger equation 
in the case of V si + V ss are very close to the masses determined by the Schrodinger equa- 
tion in the case of only V s - U a p-charmonium dissociation reaction which is endothermic 
(exothermic) in the former case is also endothermic (exothermic) in the latter case, and 
in the two cases are quite close. Since y/s dependence of the cross section for a 
quark- interchange-induced reaction is mainly determined by [30], y/s dependence of 
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the cross section for a p-charmonium reaction in the former case is similar to one in the 
latter case. Hence, we do not plot p-charmonium dissociation cross sections in the case of 
V S i + V ss , but present parametrizations of these cross sections. The cross section for the 
exothermic reaction qq(S qq ; P) + cc(S C c] —P) — > qc(S q c] P') + cq(S cq ; —P') can be related 
to the endothermic reaction qc + cq — > qq + cc by the detailed balance 

-> 2 

unpol 

_ (2S q , + l)(2S cq + l)P' 

unpol /qi \ 

a m+cc^qc+cq ~ + + l) p2 QC+cq^qq+cc, W 

where S q5 and S cq are the spins of qc and cq, respectively. It is then correct to choose the 
following parametrization for the exothermic reaction 



P 2 



+a 2 



C2 



exp 



ci | I 

'■■2 I 1 



&1 



(32) 



Parameter values are listed in Tables 4-6. 



6. Procedure 



The curves shown in Figs. 2-16 and 18-23 correspond to the zero temperature and 
the five nonzero temperatures T\ = 0.65T C , T 2 = 0.75T C , T 3 = 0.85T C , T 4 = 0.9T C and 
T 5 = 0.95T C . We now present a procedure on how to obtain the unpolarized cross section 
for qq + cc — > qc + cq for 0.65T C < T < T c from the curves. 

First, we state the procedure while a reaction at T { and T i+1 is endothermic. We 
denote by ^/s^ the square root of the Mandelstam variable corresponding to the peak 



so is the difference of ^/s^ with respect to the threshold 



cross section. d 



energy y/so~. Let be the square root of the Mandelstam variable at which the cross 
section is 1/100 of the peak cross section, and ^fsl > v /s~^ > ^/so. d , v /s^, ^/s ', and ^/Sp" 
at Tj are indicated by doi, ^J~s7i, y^oT? an d i/Spi, respectively. d 0i and can be found 
in Tables 1-6. v /s"o7 can be obtained from the mass parametrizations in Section 2, and 
V^pi = V / ^o7+ doi- For Tj < T < T i+1 (i=l, 2, 3 or 4) we take the linear interpolation 
between the two peaks at Tj and T i+ i to estimate d and ^/s^ of the cross section at T, 

(33) 



doi+i - d 0i 
d = — — (i - Li) + doi, 



Ti 



i+l 
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-(T - Ti) + y/s~. 



(34) 

v/so is the sum of final-meson masses, and parametrizations of the masses have been given 
by Eqs. (7) and (8). The square root of the Mandelstam variable corresponding to the 
peak cross section is ^/s^ = ^/sq + d . We define a ratio 



c 



if y/so < y/s < 

if ^ < y/s < y/T z . 



(35) 



The ratio ( is on the closed interval [0,1], and corresponds to y/s and T. In the cross 
section curve at T, we can find a point (y/sl, a^ npol (y/si, %)) which gives the same ratio, 



c 



if \/soi < \fs~i < 



3 pi 



(36) 



if yj s pi < y/Si < yj S z j, 

and in the cross section curve at T i+1 , we find a point (/s i+ i, a^^y/Sj+i, T i+1 )) which 



also gives the ratio 



c 



fSi+1— y/~Spi+i 



(37) 



The regions y/s^ < y/s~ < ^/s~ and y'Soi+i < v^+i - \/^pi+i (y/s& < < v 7 ^ and 
V^+T < V^+i - v^+i) correspond to y/s^ < y/s < v /s^ ( v /s^ < v/s < y^). In terms 
of the ratio y^ and y^i+i are expressed as 



/^pi + C(y / ^' 



if v 7 ^ < < 



fS~) if y/fy < yj! < y/s~ z , 



y/ s pi+l + C(v /S 0i+1 ~ -y/Spj+i) if y/Sq < \/s < 



(38) 



(39) 



y/Spi + l + ({y/S zi + 1 - y/S pi+ l) if V /S^ < < y/S^. 

The linear interpolation between the two points provides the unpolarized cross section at 
y/s for Ti < T < T i+1 , 



a 



unpol 



unpol / /- rp \ unpol / i— rp \ 

1+1 ( ^ r ' +i) - g « (^) (T _ T .) + a r po \y/r t , T) a y/r Q < y/S < y/sl 



Ti+i-Ti 



li yfs> y/s~ z . 



(40) 
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where a^ npo and cr™i° are shown in Figs. 2-16 and 18-23, and are given by Eq. (30) 
together with the parameters listed in Tables 1-6. For 0.95T C < T < T c Eqs. (33)-(40) 
still apply to obtaining the unpolarized cross section so long as Tj = T 4 and T i+ i = T 5 are 
set. 

Second, we state the procedure while a reaction at Tj and T i+1 is exothermic. Eqs. 
(33)-(40) are suited to endothermic reactions of which each has the zero cross section at the 
threshold energy or at the infinite total energy of initial mesons and has one maximum 
cross section in its y/s dependence. The cross section for the exothermic reaction is 
infinite at the threshold energy. But the quantity enclosed by the braces in Eq. (32) has 
the general y/s dependence of endothermic reactions. Hence, Eqs. (33)- (39) apply to the 
quantity enclosed by the braces, doi and of the quantity are listed in Tables 4-6. y/so 
equals the sum of initial-meson masses, and parametrizations of the masses have been 
given by Eqs. (4), (5), (6), and (24). Eq. (40) now gives the unpolarized cross section for 
the exothermic reaction at T while o-" npo1 an d cr^° l are the cross sections shown in Figs. 
8 and 11-16, and are given by Eq. (32) together with the parameters listed in Tables 4-6. 
Since the infinity of the cross section at the threshold energy is intractable, do not let y/s 
equal y/so while fortran code is made. 

Third, we state the procedure while a reaction is exothermic at Tj and endothermic at 
T i+1 . Eqs. (33)-(39) apply to the endothermic reaction at T i+1 or the quantity enclosed by 
the braces in Eq. (32) at Tj. Relevant d 0i , y/s^i, d 0i+ i, and v /s^+T can be found in Tables 
4-6. v /soi and ^/soi+i equal the sum of initial- meson masses and the sum of final- meson 
masses, respectively, ^/sq at T is the sum of final-meson masses for m q5 +m cq —m qq —m C c > 
or the sum of initial-meson masses for m q5 + m cq — m qq — m c5 < 0. After estimating 
and y/Si+i, Eq. (40) is used to get the unpolarized cross section at -/s for Tj < T < T i+i 
while (r" npol (y / ii, Tj) is the cross section for the exothermic reaction and cr^°f o1 (y'Sj+i, T + i) 
for the endothermic reaction. 

7. Summary 

The central spin-independent and temperature-dependent potential has been used in 
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the Schrodinger equation to obtain the temperature-dependent masses of J/ip, ip', Xc 
D, D* , D s and D* as well as the wave functions of quark-antiquark relative motion 
inside these mesons. The experimental masses of the mesons are reproduced by the 
potential at T = 0. While temperature increases, the temperature dependence of the 
theoretical masses is: the ip' and Xc masses decrease slowly, and the J/tp mass more 
slowly; each of the D and D s masses keeps almost unchanged in a temperature region, 
and apparently falls off near T c ; the D* and D* masses first decrease slowly, and then 
apparently fall off. The prominent medium effects are: ip' and Xc are degenerate in mass 
for 0.6T C < T < T c ; very near the critical temperature J/ip, ip', and Xc become a mass 
triplet, and D and D* (D s and D*) become degenerate in mass. These are also true when 
the central spin-independent potential and the smeared spin-spin interaction are used in 
the Schrodinger equation to obtain temperature-dependent masses and quark-antiquark 
relative-motion wave functions. The particular temperature dependence of the J/ip mass 
and radius causes the peak cross sections of the endothermic nJ/ip and pJ/ip reactions 
to increase rapidly when temperature approaches the critical temperature. Even though 
■ip' and Xc are degenerate in mass, the cross sections of ip' and Xc m collisions with a light 
meson are different, which is owed to the node of the ip' wave function. The temperature 
dependence of the potential, the quark-antiquark relative-motion wave functions, and 
the meson masses leads to the temperature dependence of all dissociation cross sections 
by the three factors, the difference m q5 + m c q — m q q — m c5 , and the transition 

amplitude. Some reactions are purely endothermic while the others become exothermic 
below certain temperatures. Parametrizations of unpolarized cross sections are given for 
nJ/ip ->■ D*D or DD*, nJ/ip ->■ D*D*, n^' ->■ D*D or DD*, nip' ->■ D*D*, n Xc ->■ D*D 
or DD*, n Xc ->■ D*D*, pJ/ip ->■ DD, pJ/ip ->■ D*D or DD*, pJ/ip ->• D*D*, pip' ->■ DD, 
pip' ->■ D*D or DD*, pip' ->■ D*D*, p Xc ->■ ^^D, PXc ->■ or Z)D*, and p^c ->■ 5*^*- 

The parametrizations at the five temperatures, T/T c = 0.65, 0.75, 0.85, 0.9, 0.95, can be 
used in the procedure to yield cross sections at any temperature between 0.65T C and T c . 
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Figure 1: Meson masses as functions of T/T c . The masses of J/ip, ip', Xc, D, D*, D s and 
D* are shown by the solid, dashed, dotted, lower long dashed, upper long dashed, lower 
dot-dashed and upper dot-dashed curves, respectively. 
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Figure 3: Cross sections for irJ/ip — > D*D* at various temperatures. 
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Figure 4: Cross sections for nip' — > D*D or DD* at various temperatures. 
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Figure 5: Cross sections for nip' — > D*D* at various temperatures. 
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Figure 6: Cross sections for nx c — >• D*D or DD* at various temperatures. 
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Figure 8: Cross sections for pJ/ip — > DD at various temperatures. 
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Figure 9: Cross sections for pJ/ip — > D*D or DD* at various temperatures. 
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Figure 10: Cross sections for pJ/ip — > D*D* at various temperatures. 
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Figure 11: Cross sections for pip' — > DD at various temperatures. 
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gure 12: Cross sections for pip' — >■ D*D or DD* at various temperatures. 



36 



0.65 



0.85 



0.95 



0.9 



0.75 



T/T c =0 



1500 



1250 



1000 



750 



500 



250 



2.8 3.3 3.8 3 3.5 4 4.5 5 



s 1/2 (GeV) 



Figure 13: Cross sections for pip' — > D*D* at various temperatures. 
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Figure 14: Cross sections for p\ c DD at various temperatures. 
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Figure 15: Cross sections for px c —> D*D or DD* at various temperatures. 
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Figure 16: Cross sections for px c —> D*D* at various temperatures. 
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Figure 17: n, K and 77 masses as functions of T/T c . 
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Figure 18: Cross sections for nJ/ip — > D*D or DD* at various temperatures. 
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Figure 19: Cross sections for nJ/ip — > D*D* at various temperatures. 
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Figure 20: Cross sections for mp' — > D*D or DD* at various temperatures. 
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Figure 21: Cross sections for nip' — > D*D* at various temperatures. 
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Figure 22: Cross sections for irx c — > D*D or DD* at various temperatures. 
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Figure 23: Cross sections for TiXc — >• D*D* at various temperatures. 
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Table 1: Values of the parameters in Eq. (30) for the irJ/ip dissociation. a\ and a 2 are in 
units of mb; b±, b 2 , do, and are in units of GeV; c\ and C2 are dimensionless. 



Reactions 


T/T c 


ai 


61 


Cl 


a 2 


b 2 


c 2 


do 




ixJli) ->■ D*D 





0.16 


0.03 


0.53 


0.05 


0.08 


2.42 


0.04 


4.556 




0.65 


0.2 


0.021 


0.51 


0.1 


0.069 


1.8 


0.03 


4.279 


or DD* 


0.75 


0.27 


0.025 


0.47 


0.1 


0.044 


0.87 


0.03 


4.089 




0.85 


0.21 


0.021 


0.5 


0.09 


0.073 


1.87 


0.03 


3.847 




0.9 


0.18 


0.015 


0.51 


0.08 


0.069 


1.75 


0.02 


3.649 




0.95 


1.05 


0.006 


0.49 


1 


0.025 


1.27 


0.01 


3.291 







0.026 


0.26 


4.54 


0.006 


0.35 


0.59 


0.25 


5.072 




0.65 


0.022 


0.2 


2.14 


0.005 


0.14 


0.56 


0.23 


4.864 




0.75 


0.026 


0.2 


3.29 


0.003 


0.52 


0.2 


0.22 


4.666 




0.85 


0.01 


0.016 


0.5 


0.006 


0.27 


5.09 


0.02 


4.287 




0.9 


0.006 


0.039 


1 


0.003 


0.005 


0.45 


0.03 


4.066 




0.95 


1.61 


0.007 


0.54 


0.45 


0.026 


2.78 


0.01 


3.187 
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Table 2: The same as Table 1 except for nip'. 



Reactions 


T/T c 


ai 




Cl 


a 2 


b 2 


c 2 


d 




TTlp' ->■ D*D 





1.78 


0.02 


0.5 


1.15 


0.06 


1.69 


0.035 


4.567 




0.65 


2.89 


0.01 


0.68 


1.37 


0.01 


0.34 


0.01 


4.021 


or DD* 


0.75 


3.18 


0.01 


1.11 


0.61 


0.01 


0.01 


0.01 


3.871 




0.85 


1.46 


0.005 


0.49 


0.89 


0.009 


0.55 


0.01 


3.606 




0.9 


1.38 


0.01 


5.59 


0.68 


0.03 


0.01 


0.01 


3.431 




0.95 


5.04 


0.002 


0.32 


3.34 


0.01 


6.61 


0.01 


3.215 







1.6 


0.025 


0.5 


0.3 


0.042 


1.66 


0.03 


4.78 




0.65 


1.6 


0.005 


0.51 


0.78 


0.018 


1.64 


0.01 


4.103 




0.75 


1.3 


0.005 


0.49 


0.37 


0.014 


1.42 


0.01 


3.916 




0.85 


0.13 


0.008 


2.36 


0.08 


0.003 


0.12 


0.01 


3.64 




0.9 


0.35 


0.005 


2.23 


0.36 


0.002 


0.28 


0.01 


3.349 




0.95 


6.4 


0.009 


5.08 


5 


0.001 


0.15 


0.01 


3.176 


Table 3: The same as Table 1 except for 7rx c - 


Reactions 


T/T c 


a\ 


bi 


Cl 


a 2 


b 2 


c 2 


d 




nXc ^ d*D 





0.95 


0.13 


3.17 


0.18 


0.07 


0.84 


0.15 


4.577 




0.65 


1.13 


0.07 


1.78 


0.5 


0.04 


0.71 


0.06 


4.135 


or DD* 


0.75 


1.39 


0.04 


1.32 


0.32 


0.05 


0.46 


0.04 


3.955 




0.85 


1.04 


0.02 


1.04 


0.16 


0.01 


0.25 


0.02 


3.636 




0.9 


1.14 


0.01 


0.83 


0.4 


0.03 


3.67 


0.015 


3.438 




0.95 


4.51 


0.017 


2.03 


1.32 


0.01 


1.05 


0.02 


3.167 







0.29 


0.08 


2.53 


0.08 


0.05 


1.12 


0.08 


5.069 




0.65 


0.39 


0.04 


1.45 


0.1 


0.05 


3.86 


0.045 


4.58 




0.75 


0.45 


0.03 


1.6 


0.11 


0.05 


1.17 


0.03 


4.293 




0.85 


0.08 


0.016 


1.48 


0.03 


0.028 


1.65 


0.02 


3.724 




0.9 


0.17 


0.003 


2.22 


0.25 


0.012 


1.98 


0.01 


3.405 




0.95 


6.36 


0.003 


2.43 


7.48 


0.015 


2.16 


0.015 


3.159 
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Table 4: Values of the parameters in Eqs. (30) and (32) for the pJ/ip dissociation. a\ and 
a 2 are in units of mb; b±, b 2 , d , and are in units of GeV; C\ and c 2 are dimensionless. 



Reactions 


T/T c 






Cl 


a 2 


b 2 


c 2 


d 




pj/i) -> DD 





0.073 


0.04 


0.51 


0.049 


0.45 


5.01 


0.05 


4.747 




0.65 


0.24 


0.008 


0.47 


0.16 


0.03 


0.79 


0.01 


4.29 




0.75 


0.22 


0.01 


0.41 


0.11 


0.01 


0.54 


0.007 


3.676 




0.85 


0.16 


0.007 


0.52 


0.049 


0.02 


0.45 


0.007 


3.887 




0.9 


0.25 


0.005 


0.37 


0.15 


0.01 


1.99 


0.01 


3.423 




0.95 


0.83 


0.005 


0.5 


0.7 


0.02 


1.57 


0.01 


3.179 


pj/ijj ^D*D 





0.71 


0.03 


0.39 


0.56 


0.42 


6.42 


0.03 


5.011 




0.65 


0.28 


0.013 


0.5 


0.067 


0.3 


2.63 


0.02 


4.648 


or DD* 


0.75 


0.15 


0.012 


0.47 


0.044 


0.27 


3.06 


0.01 


4.541 




0.85 


0.091 


0.009 


0.44 


0.045 


0.18 


3.92 


0.01 


4.196 




0.9 


0.28 


0.005 


0.38 


0.15 


0.006 


1.49 


0.01 


3.6 




0.95 


2.03 


0.007 


0.53 


0.68 


0.024 


2.18 


0.01 


3.181 


pj/ip ->■ D*D* 





0.63 


0.05 


0.47 


0.43 


0.36 


5.79 


0.04 


5.043 




0.65 


0.086 


0.06 


0.46 


0.043 


0.28 


4.05 


0.04 


4.901 




0.75 


0.043 


0.03 


0.46 


0.033 


0.24 


2.39 


0.03 


4.765 




0.85 


0.0061 


0.012 


0.46 


0.0031 


0.14 


1.6 


0.01 


4.291 




0.9 


0.27 


0.007 


0.84 


0.091 


0.002 


0.18 


0.01 


3.607 




0.95 


5.51 


0.005 


0.51 


4.47 


0.02 


1.63 


0.01 


3.189 
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Table 5: The same as Table 4 except for pip' . 



J. \AZ> CI O 


TIT 

J- / J- c 




"i 


c i 




bo 


Co 




V * z 


m i,' _i. nn 


n 

w 


02 


1 3 


72 


02 


05 

w . wo 


45 

W .tiO 


1 

W. -L 


5 355 

o . ooo 




n 65 

W .WO 


001 9 


04 


53 


0012 

W .WW ±4j 


25 

W . 4j o 


2 74 

£s. I T: 


03 

w . wo 


4 848 






D 0099 


fl 099 


n 54 


o nnn7Q 

W . WWW 1 C7 


n 97 

W.^ I 


4 83 


n D3 

w. WO 


4 37 

TI.O I 




8 . 


W. W^±c/ 


n m 


n .9 


D 047 

W . U4: / 


W. WZi J- 


9 7fi 


n m 

W . W -L 


3 .7 . 

O. O i o 




n 9 

w • u 


O . J- 


006 

W .WWW 


9 


5 57 

O . O i 


002 


38 


01 

W . W X 


3 373 

O . O 1 o 




0.95 


1.07 


0.005 


0.62 


0.26 


0.01 


0.1 


0.01 


3.171 




n 

W 


n 083 


n D5 

W . WO 


n 5 

W.O 


fl D45 


n 93 


1 8 

J. .O 


n 1 

W. J. 


5 354 

O.OOt: 




o 65 

W .WO 


017 


03 


53 

W . OO 


0043 

W . WWtiO 


3 

w . o 


8 68 

O . WO 


03 

w . wo 


4 818 


or DD* 


f) 75 

W . 1 O 


f) 097 


fl DD8 


o 46 

W.TlW 


fl 091 

W . WZj J. 


n 04 


n 74 

W. I t: 


fl m 

W. W X 


4 44 

T: .tit: 




o 85 

W . OO 


O . W 


0008 


49 


4 54 


0046 

W . WWt:W 


1 27 


005 

w . w wo 


3 568 

o . o WO 




0.9 


2.09 


0.006 


1.8 


2.09 


0.003 


0.17 


0.01 


3.385 




0.95 


0.99 


0.006 


0.77 


0.62 


0.007 


0.23 


0.01 


3.173 


p^' ->■ D*D* 





0.2 


0.03 


0.51 


0.14 


0.19 


1.37 


0.05 


5.337 




0.65 


0.78 


0.05 


3.08 


0.04 


0.0013 


0.68 


0.05 


4.535 




0.75 


10.86 


0.0009 


0.44 


13.39 


0.005 


1.05 


0.005 


3.832 




0.85 


4.79 


0.005 


0.82 


6.23 


0.002 


0.48 


0.01 


3.597 




0.9 


2.07 


0.003 


0.49 


0.72 


0.012 


0.78 


0.004 


3.394 




0.95 


2.76 


0.005 


0.47 


1.08 


0.015 


0.8 


0.006 


3.185 
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Table 6: The same as Table 4 except for px c - 



Reactions 



T/T c 



ai 



Cl 



a 2 



c 2 



do 



PXc DD 




0.65 
0.75 
0.85 
0.9 
0.95 



0.029 
0.007 
0.001 
0.03 
2.98 
0.52 



0.07 
0.21 
0.007 
0.003 
0.007 
0.009 



0.5 
3.1 
0.51 
0.63 
1.42 
1.54 



0.0073 
0.0003 
0.0008 
0.0062 
1.94 
0.32 



0.21 
0.005 
0.24 
0.01 
0.017 
0.019 



1.24 
0.51 
3.59 
0.2 
2.98 
3.15 



0.1 
0.2 
0.01 
0.005 
0.01 
0.01 



PXc -)• D*D 



or DD* 




0.65 
0.75 
0.85 
0.9 
0.95 



0.1 
0.0029 
0.0034 
1.75 
1.58 
0.34 



0.16 
0.008 
0.007 
0.003 
0.01 
0.014 



1.16 
0.44 
0.5 
0.5 
1.4 
1.33 



0.045 
0.018 
0.0044 
1.41 
0.81 
0.11 



0.02 
0.23 
0.2 
0.011 
0.015 
0.013 



0.49 
2.7 
3.09 
1.63 
2.03 
4.84 



0.15 
0.25 
0.15 
0.005 
0.01 
0.01 



PXc -> D*D* 




0.65 
0.75 
0.85 
0.9 
0.95 



0.31 
1.91 
8.57 
2.21 
0.29 
0.11 



0.27 
0.005 
0.006 
0.007 
0.012 
0.013 



3.43 
0.52 
0.58 
1.58 
1.66 
1.99 



0.29 
1.87 
8.57 
1.96 
0.13 
0.041 



0.04 
0.024 
0.018 
0.018 
0.013 
0.025 



0.54 
1.55 
1.82 
3.48 
1.02 
0.76 



0.25 
0.01 
0.01 
0.01 
0.01 
0.01 
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